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ABSTRACT 

The high state of black hole candidates is characterized by a quasi-thermal emission component 
at kT ~ 1 keV. In addition, this state tends to have very low variability which indicates that it is 
relatively stable, at least on short time scales. Most models of the high state imply that the bulk of 
the emission comes from an optically thick accretion disc; therefore, this state may be an excellent 

■^j- laboratory for testing our ideas about the physics of accretion discs. In this work we consider 

the implications of assuming that accretion disc viscosity arises from some form of turbulence. 
Specifically, we consider the simple case of three dimensional hydrodynamic turbulence. It is found 
that the coupling of such turbulence to acoustic modes in the disc can alter the disc emission. We 
calculate the amplitude and frequencies of this modulation, and we express our results in terms of 
the X-ray power spectral density. We compare our calculations with observations of the black hole 
candidate GS 1124-683, and show that for certain parameters we can reproduce some of the high 

^vq frequency power. We then briefly explore mechanisms for producing the low frequency power, and 

note the difficulty that a single variability mechanism has in reproducing the full range of observed 
variability. In addition, we outline ways in which future spacecraft missions - such as USA and 
XTE - can further constrain our model, especially at frequencies above ~ 10 2 Hz. 
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6 1 INTRODUCTION 

H 

^ Black hole candidates (BHC) are typically identified either from estimates of their mass, via a 

binary mass function, or by analogies with other BHC where one typically compares high energy 
emission and variability properties. (For a review of black hole candidates, including a summary 
of observed spectral and temporal properties, see Tanaka & Lewin 1994.) The energy spectra of 
BHC have been historically labelled based upon observations of the soft X-ray band (~ 1 — 10 keV). 
Intense, quasi-thermal flux is referred to as the "high" state. Non-thermal flux (typically a power 
law with a photon index of ~ 1.7) in this band indicates that the BHC is in the "off" state, for 
extremely low intensity flux, or "low" state, for moderate intensity flux that in the ~ 1 — 10 keV 
band is less than the high state flux in the same band. [It is possible, however, for the bolometric 
luminosity to be greater in the low state than in the high state. See the comments by Tanaka & 
Lewin (1994) and Nowak (1994b).] The high state tends to have little variability, with a root mean 
square (RMS) variability of a few percent, whereas the low state tends to have an RMS variability 
of several tens of percent. (We will make explicit what we mean by RMS variability in §4.) 

Many models have been proposed to explain the observed spectra. Ultimately, nearly all of 
them involve the compact object accreting matter via an accretion disc. For the case of the non- 
thermal low state, this disc has been modelled as being optically thin (cf. Melia & Mesra 1993, 
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Luo & Liang 1994), as a cold reflector of a central power law source (cf. Done et al. 1992), and as 
shrouded in a hot, optically thick, Comptonizing cloud (cf. Sunyaev & Titarchuk 1980). On the 
other hand, the most commonly accepted model for the quasi-thermal high state is optically thick 
emission from a relatively stable and quiescent disc (cf. Shakura & Sunyaev 1973). Though the 
details of the high state solutions differ, the one common feature that they have is the emission 
domination of the disc itself. That is, there is no need to invoke clouds, winds or other external 
sources of emission. Recent models have been very succesful at explaining luminosity variations in 
this state as arising from accretion rate variations in a disc of constant inner radius, lending further 
support to the accretion disc picture (cf. Ebisawa et al. 1993, Ebisawa 1994). In addition, the 
fact that the RMS variability of the high state is very low suggests that the disc is in a relatively 
steady state, at least on short time scales. These properties of the high state - predominantly 
disc-dominated and steady accretion - are what make it the ideal laboratory for testing our notions 
about the physics of accretion discs. 

The purpose of this work is to explore how one aspect of accretion disc physics, namely disc 
turbulence, may manifest itself in the observed properties of BHC high states. Specifically, we will 
explore the ways in which turbulence affects the observed variability of the high state. In the next 
section, we will outline our basic assumptions about the nature of accretion disc turbulence. For 
our purposes, we shall only consider isotropic hydrodynamic turbulence. (We will explore the more 
realistic cases of anisotropic and/or magnetohydrodynamic turbulence in future works.) In §3 we 
will consider two ways in which turbulence may affect the observed disc emission, namely, advection 
of photons from the hot disc midplane and direct modulation of the emission rate. The latter will 
be seen to be the more important of these two. This manifestation of turbulence is essentially the 
coupling of acoustic modes to the turbulent eddies. In §4 we consider how the emission modulation 
is revealed in the observed X-ray power spectral density of BHC high states. We compare our 
predictions to high state observations of GS 1124-683 that cover the frequecy range 10~ 2 — 10 2 Hz. 
We will see that for certain parameter ranges, acoustic mode coupling to turbulence can account 
for the high frequency power. In §5 we briefly explore whether weak fluctuations on viscous or 
thermal time scales can explain the observed low frequency power. Finally, in §6 we summarize our 
results and discuss both future theoretical and observational prospects. 

2 BASICS OF THE DISC TURBULENCE 

The single greatest uncertainty in modelling accretion discs is the viscosity prescription. This 
uncertainty directly affects our understanding of mass and angular momentum transport in the 
disc, but more importantly for our purposes it affects our understanding of the energy generation 
(and hence structure) within the disc. Most models require a phenomenological description of 
the radial energy generation profile and disc structure; however, here we require a model of this 
structure in the vertical direction (i.e., perpendicular to the disc midplane) as well. In this work 
we shall explore the consequences of assuming that viscosity is ultimately due to accretion disc 
turbulence. 

One of the original motivations for the a-disc model of Shakura & Sunyaev (1973) was con- 
sideration of (locally) homogeneous and isotropic hydrodynamic turbulence. Applications and 
modifications of the a-disc model have grown beyond this initial motivation; however, in this work 
we shall return to the viewpoint that viscosity is directly related to the properties of turbulent 
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eddies. Even within this framework, recent works have considered models more sophisticated than 
homogeneous and isotropic hydrodynamic eddies. Some researchers have carefully considered the 
role of rotation and have allowed for anisotropic (i.e. two-dimensional) eddies (cf. Narayan et al. 
1994, Kato & Yoshizawa 1993, Dubrulle 1992). One of the more exciting developments has been 
the study of the magnetic shearing instability (MSI), which leads to three-dimensional magneto- 
hydrodynamic turbulence (Balbus & Hawley 1991, Hawley & Balbus 1991, Balbus et al. 1994). 
Similar to classical isotropic hydrodynamic turbulence, the MSI turbulence is quasi-incompressible. 
However, its relation to energy generation in a disc, and its coupling to the global dynamics of the 
accretion disc, are as of yet unknown. In this work we will only consider isotropic, hydrodynamic 
turbulence as embodied in the a-disc model. We choose this model as being the most analyti- 
cally tractable and as being able to provide a guide for using more sophisticated models (namely 
anisotropic eddies and MSI turbulence) in future work. For the basic physics of this model we will 
follow the work of Shakura & Sunyaev (1973) and the discussion of turbulence found in Landau & 
Lifshitz (1987). 

At a given radius in the disc, let the largest turbulent eddy size be L, and have an eddy velocity 
ul, such that 

u L ~ M L c s , L ~ Sh , (2.1) 

where Ml is the Mach number of the turbulence (assumed Ss 1), c s ~ hQ is the disc speed of sound, 
h is the half-thickness of the disc, Q is the Keplerian rotation velocity, and S is a dimensionless 
scaling factor. In the vertical direction the eddy size is limited by the disc thickness. In the radial 
direction, constraining the disc velocity difference across the eddy to being less than the speed of 
sound also limits the eddy size to Ss h. The turbulent viscosity is found by taking the simplest 
combination of the above quantities that has the dimensions of a kinematic viscosity. We therefore 
have for the turbulent viscosity 

v t ~ u L L ~ M L 8h 2 {l ~ ah 2 Q, . (2.2) 

In the a-disc model the viscous azimuthal stress is given by ~ v t p£l ~ aP, where P ~ pc\ is the 
pressure. This gives us the above identification a ~ Ml§- 

Assuming that the standard form for the viscosity tensor holds for our Keplerian disc, the rate 
of viscous energy dissipation per unit mass is given by 

5fT 2 



e ~ v t (j-j^) ~ u L LW . (2.3) 

However, from the standard Kolmogorov turbulence scaling laws (cf. Landau & Lifshitz 1987) we 
have 

u L ~ (e£) 1/3 ~ (u L L 2 Q 2 ) 1/3 =^ ul ~ LQ =^ 8 ~ Ml ■ (2.4) 
We can then identify 

a ~ M 2 L . (2.5) 

Since a is a direct measure of the Mach number of the turbulence, it is directly related to the 
variations in the energy generation rate, as we will show in subsequent sections. 
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Using a as a measure of the turbulent Mach number of the largest eddies, we find the turbulent 
velocity, ul, eddy size, L, and eddy lifetime or rollover time, t^, of the largest scale eddies to be 
given by 

ul ~ \fah£l , L ~ \fah , ~ L/ul ~ • (2-6) 

On smaller scales, I, we use the standard Kolmogorov scaling laws to find 

1/3 //>\2/3 



aV*m (£) , h ~ l/ Ul ~ a-'/'Q- 1 (£) , (2.7) 



where I Ss ^/ah. Note that for the smaller scales, the turbulent Mach number Mi ~ M 2 / 3 (i/h) 1 ^ 3 . 
The pressure variation associated with the largest eddy is given by 

Api, ~ pu\ ~ aph 2 tl 2 , (2.8) 

and therefore on smaller scales, I, we have 

/ 1 \ 2 / 3 

Ap £ ~ ~ a 2/3 ph 2 Sl 2 l-\ . (2.9) 

We have used the scaling laws for incompressible turbulence; however, we expect weak density 
fluctuations to be associated with these eddies. The density fluctuations should be given by Ap ~ 
cJ 2 Ap, therefore we have 

Ap L ~ M\p , A Pl ~ M\p • (2.10) 

We can consider these density fluctuations as being due to acoustic modes coupled to the turbulence 
(cf. Goldreich & Kumar 1988). Essentially, eddies of frequency fi = tj 1 will couple to acoustic 
modes of the same frequency or lower. In this work, we shall calculate the coupling to azimuthally 
symmetric (m = 0) modes, as they are observationally most relevant. [Modes with multiple nodes in 
the azimuthal direction (m > 1) will have the azimuthally integrated flux reduced to some degree.] 
We shall use equation (2.10) as an estimate of the amplitude of this coupling for the calculations 
of §4. 

3 INFLUENCE OF TURBULENCE ON THE DISC SPECTRUM 

In this section we shall consider ways in which the accretion disc turbulence might manifest itself 
in the observed X-ray variability of the source. We consider two possible processes: advection 
and direct modulation of emission. If the typical photon radiative diffusion time from the disc 
midplane is much longer than a typical eddy rollover time, there exists the possibility that eddies 
will advect photons from the hot disc midplane to the cooler disc surface. (For a Shakura-Sunyaev 
a-disc model, the midplane temperature is typically a factor of ~ 10 greater than the disc surface 
temperature, though this depends upon such parameters as a and the accretion rate; Shakura & 
Sunyaev 1973, Shapiro & Teukolsky 1983.) One must also check to see that the photon absorption 
time is longer than the eddy rollover time, else the photon will be absorbed and reemitted with 
properties of the local temperature. 

For the case of emission modulation, if one assumes that the energy generation rate per unit vol- 
ume scales as the local density, then the emission can be modulated by the weak density fluctuations 
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associated with the turbulence. In this case, the photon diffusion time (though not necessarily the 
photon absorption time) must be longer than the eddy lifetime in order for the over /under- density 
to be manifested in emission variability. We shall explore each of these possibilities separately 
below. 

a) Advection 

In order for advection to be effective at "dredging" up photons from the hot midplane of the disc, 
the photons must undergo many scatterings during the eddy rollover time. We assume that these 
scatterings will occur as the result of electron scattering opacity, n es = 0.4 cm 2 g _1 . The photon 
diffusion time due to scattering across an eddy length, I, must be much greater than the rollover 
time of the eddy, t^. (Equivalently, the photon diffusion length over a time ti must be much less 
than I.) We shall ignore the effects of absorption in the following, which would serve to make our 
constraints even more severe. 

The number of scatterings, Ni, over a diffusion length, I, is that of a random walk (cf. Rybicki 
& Lightman 1979), namely 

n\' 2 ~ l/\ es ~ K esP l = T e t s , (3.1) 

where A es is the electron scattering mean free path and r| s is defined to be the electron scattering 
optical depth over length I. The total path length travelled by the photon is given by 

N t \ es ~ N\' 2 l ~ rfl = T e h s h (J^j , (3.2) 

where t^ s is the electron scattering optical depth over the disc half-thickness h. The diffusion time 
over length I is therefore 



where c is the speed of light. For advective transport to be significant, this diffusion time must be 
much greater than the eddy rollover time: 

Accounting for the maximum eddy size as given in equation (2.6), the above implies that 

^ J> (£) » a" 1 / 4 KT 3/i (^)" 3/4 • (3.5) 

In a previous work (Nowak & Wagoner 1992) we have noted that for a Shakura-Sunyaev a-disc, 
c s /c ~ /3 _1 >C, and that £ = 2ph ~ /3 2 (a K es £) _1 , where /3 is the ratio of radiation pressure to 
total pressure and C is the ratio of the disc luminosity to the Eddington luminosity of the central 
compact object. Using these estimates, we find that r^ s c s /c ~ /3/a, and therefore for photon 
advection to be effective in advecting photons from the midplane to the surface (giving the greatest 
temperature contrast) we require 

a % » f3~ 3/i V^ ■ (3.6) 
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Since /3 ^ 1, this effectively excludes photon advection as an important process for modulating the 
disc output. 

b ) Emission 

The unimportance of advection can also be an advantage. The fact that a photon will escape the 
turbulent eddy in which it was created in less than an eddy lifetime means that photons will retain 
knowledge of the eddy's over /under- density. That is, modulation of emission will not be smeared 
out by diffusion over many eddy lifetimes. The problem of emission modulation breaks up into two 
parts - the atmosphere above effective optical depth r* ~ ^T abs (r abs + r es ) ~ 1 (which produces 
an optically thin spectrum and has the potential to effect the shape of the emergent spectrum via 
electron scattering); and the atmosphere below an effective optical depth r* ~ 1 (where the bulk 
of the energy generation occurs). For both parts of this problem, modulation from the small scale 
(high frequency) eddies at large depths will be smeared out because photons diffuse through the 
remaining atmosphere above the eddy (typically a greater length than the eddy size itself) in times 
longer than the eddy lifetime. 

First let us consider emission from depth r* ^ 1. The diffusion time from a vertical depth, z, 
is given by 



j.D es 



h l z\ 2 , 

-Ah) ■ (3J) 



In order that any modulation not be suppressed, we require Ss ti (cf. equation[3.4]). This 
implies that for modulation due to an eddy of size I we are restricted to depths of 

(i) , „-./. {<) - . (3 . 8) 

With our previous estimates of the disc parameters, this becomes 

(!) s ^""-'"{ff" • < 3 - 9 > 

(Note that the above depth is always greater than an eddy size, in agreement with the results of 
the previous section.) As an example, for typical accretion disc parameters (a ~ 0.1, M ~ M@) 
effective optical depth r* ~ 1 occurs at z/h ~ 0.01 — 0.02, which implies that in order for an eddy 
to modulate the emission, it must have 

y J> 10" 4 -10- 5 . (3.10) 
h 

This covers a fairly broad range of eddy sizes and frequencies; however, being that this is an effective 
optical depth r* ^ 1, the modulation will not affect the bulk of the emitted energy which mainly 
originates at optical depths r* <; 1. 

We now consider the more important mechanism of diffusion from the energy generating part 
of the atmosphere at r* > 1. The diffusion time outward from the midplane of the disc is given by: 

r h t max ~ r h r h - 1—^-1 ~ , (3-11) 
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which yields 



where t^ ax is the diffusion time over a distance Az maa . equal to the thickness of effective optical 
depth unity, and is the total effective optical depth of the atmosphere. Implicit in the above is 
the condition 1 < < r^ s . We require the diffusion time from the midplane to be less than an 
eddy lifetime, implying 

-- 1 ' 3 {%)" * fl • < 312 > 

I * (f C j) 3 " ■ (3.13) 

Eddies of this size are able to modulate the emission on their rollover time scale even at the greatest 
depths. 

We can evaluate the above inequality for the X-ray producing regions of Shakura-Sunyaev a- 
discs. These discs can be broken up into either radiation pressure dominated regions or gas pressure 
dominated regions. Utilizing the disc solutions of Nowak (1992) (based upon those of Shakura & 
Sunyaev 1973), the above inequality becomes 

Y £ v^a (1 " 1000) C? , (3.14) 
h 

in the radiation pressure dominated regime, and 

L - J> (10- 3 - 5) £ 21 / 20 , (3.15) 

h 

in the gas pressure dominated regime. That is, as long as C Sa 0.1, at least the largest eddies 
(with l/h ~ y^a) will be able to modulate the emission. Over a wide range of radii and shallower 
depths, smaller (and hence higher frequency) eddies will also be able to modulate the emission of 
an accretion disc. In the next section, we explore the amplitude and form of this modulation for 
accretion discs about galactic black hole candidates. 

4 RELATIONSHIP TO THE POWER SPECTRAL DENSITY 

We have seen above that accretion disc turbulence will mainly affect the observed flux by coupling 
to acoustic modes which in turn alter the emission from various radii in the disc. The question 
now becomes one of how these acoustic modes manifest themselves in the observed variability as 
characterized by the normalized X-ray power spectral density (PSD). Given an observed photon 
number count rate, x, the PSD is formed from the square of the Fourier transform of a;. Integrating 
the PSD over all Fourier frequencies yields a quantity that is proportional to (x 2 ) , where the brackets 
indicate an average over the photon count rate time-series. For the particular PSD normalization 
that we shall use, integrating the PSD over / > yields: ((a: 2 ) — {x) 2 )/{x) 2 (the square root of this 
quantity is the RMS variability, cf. Miyamoto et al. 1992). Note that with this normalization the 
PSD is independent, to within noise, of source luminosity and distance. 

We assume that the photon intensity from the disc is directly proportional to the disc's local 
energy generation rate. That is, we shall ignore variations in the shape of the emergent photon 
spectra. (This assumption is equivalent to saying that locally AlnX >• AlnT, that is the fractional 
change in the luminosity from a given radius is greater than the fractional change in the temperature. 
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For a blackbody spectrum, this assumption is essentially the same as saying AlnT 4 ~ AlnT 3 , 
which we consider to be accurate enough for our crude calculations.) Furthermore, we assume that 
the disc's energy generation rate scales as the local density to a power of 0(1). This assumption is 
implicit in many a-disc models (cf. Shakura & Sunyaev 1973), though some models assume that 
the energy generation rate is proportional to the local pressure (cf. Laor & Netzer 1989). The 
two assumptions are equivalent in the gas pressure dominated regions of the disc, but not in the 
radiation pressure dominated regions. For our purposes, we take the variations in the photon flux 
to be directly proportional to the variations in the density. 

The mean square photon count rate, (a: 2 ), is composed of two parts: the mean (DC) intensity 
squared, plus an incoherent sum of oscillating acoustic sources. Thus we have: 

' (4-1) 



Jv Jf \ frequency per volume J J 



where (Aa; 2 ,) is the mean square photon count rate due to a single mode and where the integrations 
are over the volume of the accretion disc and over the frequencies of the acoustic modes. For each 
individual mode we can write 



where Cp = photons/sec/volume, or CpAV ~ Aa;. The constant, C, is related to the disc photon 
flux T via 

dx 

where r is the radial coordinate of the disc and h is the disc half-thickness. We therefore have for 
the mean square photon count rate 

,^ = ^ + / / v/ / /(||)^^)<l^] > , (4.4) 

where f~ 2 d 2 N/dVdf is the acoustic mode phase space density, and v(f) is the acoustic mode 
volume. Note the similarity of the above to the previous definition of the integrated PSD. We 
therefore identify as the normalized PSD 





dr2 „ inm^ii^ 



I 

W) = x —rf v x ' ■ (4-5) 

( / dr 2-Kr F(r) j 

In the above definition of the normalized PSD, the unperturbed accretion disc model provides 
J-(r) and h, the model of turbulence provides Ap/p, and the nature of the acoustic modes provides 
d 2 N/dVdf and v(f). We shall only study the case of axially symmetric modes as these will be be 
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the most readily observable (being of uniform brightness around a ring, there are no cancellations 
in the observed flux modulation). The dispersion relation of acoustic modes in a rotating disc is 
given by 

a; 2 « c\k 2 + k 2 , (4.6) 

where w is the mode frequency, c s is the speed of sound, k is the mode wave number, and k is 
the epicyclic frequency (k 2 = 4f2 2 + rdQ 2 /dr, where Q is the Keplerian frequency, cf. Nowak & 
Wagoner 1991). In terms of the mode frequency, /, given in Hz, this can be written as 



/« ^ + A,-) . (4.7) 

With this dispersion relation, the mode volume becomes 

•-^-""(tf (t)'- '"*'(£)"' (£)' ' (4 ' 8) 

In the above, fi, = Q is the local Keplerian frequency, i.e. the frequency of the largest turbulent 
eddies, and fi is the rotation frequency, in Hz, of the disc inner edge. Taking the frequency 
derivative of the inverse of the mode volume, we find the phase space density to be 

<PN d ! 1 fX r X z \ f 

(4-9) 




(Note that in equations [4.8] and [4.9] we have taken a "worst case scenario" of A,, ~ X z ~ h.) 

For the model of acoustic turbulence that we have chosen, the density fluctuation, and hence 
the local photon count rate fluctuation, depends upon 

^£ ~ MJ ~ ^"(ff" ~ I . (4.0) 

Thus as a function of frequency, /, we have 

2 





(4.11) 

Below we briefly will explore other frequency dependences for the density fluctuation, but we will 
see that it makes no difference in the observable part of the PSD. 

Note that in equations (4.8-4.11) we wrote the quantities in terms of fi, / fi, the local Keplerian 
frequency divided by the Keplerian frequency of the disc inner edge. We do this because the one- 
to-one correspondance between disc radius, r, and Keplerian frequency, f^, allows us to replace the 
radial integration in equation (4.5) with a frequency integration. Performing a change of variables, 
we replace r and dr with 

-2/3 , v -5/3 , v 

and dr 4 ( ^ ) d ( ^ ) . (4.12) 
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We must keep in mind that at each radius, the local Keplerian frequency is the lowest turbulent 
frequency. That is, locally turbulence will excite acoustic modes with frequencies greater than or 
equal to the Keplerian frequency. If we consider the PSD at a given frequency, /, it will be made 
up of acoustic modes that range over radii from r — > oo to a radius r such that fl(r) = f (for 
/ < fi, otherwise r cuts-off at the disc inner edge). We can therefore write for the normalized PSD 

V(f) ~ 3 X !0- ^ £ (|) - f ° </ "" 1, V (f ) 4 (f ) . (4. 13 ) 

In the above, Af is the normalization integral 

Af = (^J drr F(r)j . (4.14) 

Furthermore, for the integrations involved in equations (4.13-4.14), r and h are now in units of 
GM/c 2 (G is the gravitational constant, M is the central object mass, and c is the speed of light) 
and T can be considered completely dimensionless, having all dimensional constants factored out 
from the top and bottom integrals. 

As an example, let us look at the case of an optically thick, gas pressure dominated disc. 
For this case, the photon flux is a function of the effective temperature, which in turn is a simple 
function of the radius. We find 




T 2 oc T! ff oc r-" 9 / 2 ~ 6" 9 / 2 [ ^ ] . (4.15) 
For the disc thickness we have 

(\ -7/io 

^ oc r- 21 / 20 , (4.16) 

where 

h ~ 0.1 C 1 ' 6 a- 1 ' 10 M- 1 ' 10 , (4.17) 

and C, as before, is the disc luminosity expressed as a fraction of the Eddington luminosity and AA 
is the central object mass in units of Mq (c/. Shakura & Sunyaev 1973). Thus we find 

i ^/BqIB/IO / f \ " 7 /io 

W) ~ 3 x 10-* /r 1 T > ( 4 - 18 ) 



Mi/" y f . 

for f < fi- The frequency on the inner edge of the disc is /; ~ 2.2 X 10 3 /.M Hz. This is well beyond 
the range of most detectors, so we choose to normalize the PSD to a frequency of 10 2 Hz. We then 
have 

for f < fi- (Note that more realistic flux laws rollover near the inner edge of the disc and hence 
produce fewer photons in the inner regions than implied by eq. [4.15]. This will tend to decrease 
Af and hence increase the amplitude of th PSD by a factor of several.) 
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The slope and amplitude of the PSD are in very good agreement with the data obtained from 
GS 1124-683 as presented in Figure 1. The data in this figure are taken from Miyamoto et al. 
(1993) and show the V(f) oc f~ 7 characteristic of black hole high states. Note that for f < fi 
the theoretical PSD slope is entirely determined by the radial dependence of the disc thickness, 
whereas for f > fi the slope is determined by the physics of the turbulent cascade. Specifically, if 
we assume h oc r@ and that for the turbulent cascade the turbulent kinetic energy per wavenumber, 
E(k), scales with wavenumber, k = I -1 , as E{k) oc k~ v , we obtain 



V(f) oc 



a 2 h ( f 



fi V fi 



-2/3/3 



OC 



fi \ /i 



I ]<1 



(4.20) 



For the Kolmogorov cascade v = 5/3, and therefore the high frequency power rolls over as /~ 5 . 

That the slope and amplitude from equation (4.19) match the observations so well is most 
likely merely fotuitous, as we have not included the detector bandpass. Equation (4.19) counts 
photons of all energies, whereas the data of Figure 1 is for the bandpass 1.2 — 15.7 keV. Including 
this bandpass will tend to exclude the low Keplerian frequency outer regions of the disc, which 
emit predominantly low energy photons. Below, we explore two accretion disc models for which we 
compute the effects of including the bandpass. 

The first model is a variant of the Shakura-Sunyaev a-disc model (Nowak 1992). The disc is 
in Keplerian rotation starting from the marginally stable orbit at r = 6 out to infinity (throughout, 
we shall take the radius, r, and disc half-thickness, h, to be in units of GM/c 2 ). The flux from a 
given radius is derived by assuming energy conservation as mass flows from one Keplerian rotating 
ring to another. The energy flux from both sides of the disc is then given by 



F(r) = 1.24 X 10 28 erg cm -2 s _1 — 1 1 r~ 3 



(4.21) 



where C and M. are defined as before, 



6 36 

1 ~- + ^2 



8 60 

1 ~- + ^2 



(4.22a) 



and the function 




(4.226) 



arises as a consequence of assuming that the viscous torque vanishes at the disc inner edge (cf. 
Shakura & Sunyaev 1973). The disc is broken into radiation pressure dominated and gas pressure 
dominated regions. (For the most part we will ignore the issue of disc instabilities, though see the 
discussion of §5.) In general, the radiation pressure dominated regions are located in the inner 
regions of the disc, and the width increases as the disc luminosity increases. In the radiation 
pressure dominated regions we take the disc half-thickness (in units of GM/c 2 ) to be 



h 



18. C 1 



(4.23) 
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while in the gas pressure dominated regions we take 



h « 1.7 Xl0" 2 ( — ) j 1 / 5 r . 21 / 20 . (4.24) 
\aM J v y 

(The transitions between the regions are easily calculated from the disc structure equations; cf. 
Shakura & Sunyaev 1973, Shapiro & Teukolsky 1983.) 

Since we are assuming that the disc is optically thick throughout, the emitted spectrum will be 
either a blackbody spectrum, or a modified blackbody spectrum (cf. Shapiro & Teukolsky 1983). 
For the blackbody spectrum, we assume that at each radius the disc emits a blackbody spectrum 
at a temperature T e f f given by 

kT eff = 10.5 keV (^) 1/4 f 1 / 4 7 1-3 / 4 I 1 / 4 . (4.25) 

Note that this temperature peaks at a value kT e ff ~ 1 keV (C/M) 1 / 4 , which is typically a lower 
value than is observed for the peak temperature of black hole candidate systems in their high state. 
In the electron scattering opacity dominated regions of the disc, we expect the emitted spectrum 
to take on a modified blackbody form, for which the flux per unit frequency is given by 

z 3 / 2 exp(-z/2) 
« [exp(z) - !]■/* • < 4 ' 26 > 

where z = hv /kT s , and T s is the modified blackbody temperature (cf. Shapiro & Teukolsky 1983). 
For the gas pressure dominated regions of the disc this temperature is given by 

kT s « 2.5xl0 3 keV ( J 8 / 9 r ~ 5 / 3 , ( 4 .27) 

while for the radiation pressure dominated regions of the disc we take 

kT s « 24keV (^h-j j 16 / 45 f.- 87 / 90 . (4. 28) 

In practice, we use a modified blackbody spectrum whenever kT e ff < 2.25 kT s . (At this transition 
point, the bolometric photon fluxes of the blackbody and modified blackbody spectra are identical.) 

Using the disc structure equations, we self-consistently calculate the transitions between radi- 
ation and gas pressure dominance, and blackbody and modifed blackbody spectra. The positions 
of these transitions, as well as the value of the various disc temperatures, are dependent upon a, 
C, and Ai. Many of the black hole candidates with well determined mass functions appear to 
have masses near M ~ 6 M@, therefore we shall set M. = 6. Furthermore, for the high, quiet 
state of black hole candidates, C ~ 0.1 (cf. Nowak 1994b); therefore, we shall explore the param- 
eters C = 0.03, 0.1, 0.3. At each of these luminosities, we calculate a PSD for the parameters 
a = 0.03, 0.1, 0.3, 1.0. We use the same bandpass as was used for the GS 1124-683 data presented 
in Figure 1, namely 1.2 — 15.7 keV. The bandpass is included in equations (4.13) and (4.14) by 
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multiplying the photon number flux, J-{r), by the fraction of photons emitted into the bandpass at 
that radius. Results for these 12 runs are presented in Figure 1. 




1Q -2 1Q -1 1Q 1Q 1 1Q 2 1Q 3 1Q -2 1Q -1 1Q 1Q 1 1Q 2 10 3 




1Q -2 1Q -1 1Q jqI 1Q 2 1Q 3 1Q -2 10 -1 jqO 1Q 2 1Q 3 

Frequency (Hz) Frequency (Hz) 



Figure 1. Comparison of theoretical acoustic mode PSD with GS 1124-683 data (circles; cf. 
Miyamoto et al. 1993). In all graphs, solid line: a = 1, dotted line: a = 0.3, short dash: a = 0.1, 
long dash: a = 0.03. The disk models were as follows. Upper left: Shakura & Sunyaev disk with 
30% of the Eddington luminosity. Upper right: Shakura & Sunyaev disk with 10% of the Eddington 
luminosity. Lower left: Shakura & Sunyaev disk with 3% of the Eddington luminosity. Lower right: 
Blackbody at each radius with T 4 proportional to the local flux and the peak temperature artificially 
set to 1 keV at 10% Eddington luminosity. For all models, the central object mass was 6 M@, and 
the bandpass was 1.2 — 15.7 keV (except for the long dotted line of the lower right figure which 
included all photons). 

We compare these results to the PSD of GS 1124-683 in its high, quiet state. For these 
observations, C ~ 0.1 (Miyamoto et al. 1993). There are several things to note here. First, 
as expected, including the bandpass drastically reduces the low frequency power. The power is 
reduced to such an extent that none of the models have as much power as the observations below 
/ ~ 10 Hz. Second, we require a fairly large a in order to fit the high frequency observations. 
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This is because larger a implies both larger acoustic mode amplitudes (cf. equation [4.11]) as 
well as higher disc temperatures that give significant flux in our bandpass (cf. equations [4.27], 
[4.28]). Setting a = 1.0 clearly gives too much power, a = 0.3 is in reasonable agreement with the 
high frequency observations, and a < 0.3 for the most part gives too little power. Unfortunately, 
the observations are below that of the expected rollover in the acoustic mode power. If the high 
frequency power is really due to disc turbulence, then high frequency power spectra - such as should 
be available with the upcoming experiments on the XTE and USA satellites - should show a sharp 
drop in the PSD above / ~ 10 2 Hz. 

Since the emitted spectrum of an accretion disc is highly uncertain, we explore one other 
accretion disc model. For this model we choose the energy flux to have the same functional depen- 
dence as in equation (4.21); however, we assume that the disc emits as a blackbody throughout 
but artificially set the peak temperature to 1 keV. We choose 1 keV as being typical of the obser- 
vations, and take T e ff oc F 1 / 4 as being the least parameter dependent approximation (there is no 
dependence upon the unknown parameter a, and by artificially setting the peak temperature to 
1 keV we have removed the dependence upon C and Ai). We calculate the PSD for this model 
both with the 1.2 — 15.7 keV bandpass and without it. For these calculations we set C = 0.1, 
M. = 6, and a = 0.3. All three parameters determine the radii at which the disc makes a transition 
between radiation pressure dominance and gas pressure dominance, and furthermore a determines 
the amplitude of the acoustic modes. The results are also presented in Figure 1. Note that the 
simulation with no bandpass agrees with the data quite well; however, adding in the bandpass cuts 
off the low frequency power. There is reasonable agreement between the simulation and the data 
for / ~ 3 — 10 2 Hz, which is a slightly broader region of agreement than the Shakura-Sunyaev disc. 
However, there is still two and a half orders of magnitude of low frequency data that cannot be 
explained by acoustic modes. 

5 VISCOUS/THERMAL EFFECTS 

As we have seen from the previous section, fluctuations on turbulent (dynamical) time scales can 
only account for the observed power between a few and several hundred Hz. Fluctuations on time 
scales ranging from ~ 10~ 2 — 3 Hz must arise from other sources. The purpose of this section is to 
show that a plausible source of this low frequency power is weak variations of the disc flux on local 
viscous/thermal time scales. We expect that such fluctuations might arise from disc instabilities (cf. 
Shakura & Sunyaev 1976, Piran 1978). The standard Shakura-Sunyaev a-disc becomes unstable to 
thermal and viscous instabilities at C ~ 0.03 for typical parameters. If C <; 0.1, nearly the entire 
X-ray emitting region of the disc is subject to viscous and thermal instabilities. Below we will 
not consider any detailed mechanisms for disc instabilities; rather, we will merely show that weak 
flux variations on characteristic instability timescales can plausibly account for the low frequency 
power. 

We assume that locally the accretion rate is only steady when averaged over times greater than 
the local viscous/thermal time scales. For times shorter than this we assume that the flux varies 
exponentially in time with some fractional maximum amplitude, _;, over a radial extent, Ar. The 
variation in the photon count rate, Aa:, emitted from radius r then becomes 

Aa: = 2tt rAr jF . (5.1) 
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We know that integrating the power spectral density over frequency yields a quantity proportional 
to the square of the above variation. Furthermore, exponential variations yield PSDs that have 
functional forms V(f) oc 1/[1 + (27r/r) 2 ], where r is the viscous time scale. We therefore adopt this 
as the functional form for the PSD at a given radius (cf. Nowak 1994a). We find the amplitude, 
A, of the PSD of this radius by setting 



Ax 2 

which yields 



A = 16tt 2 r r- 2 Ar- 2 j 2 T 2 . (5.3) 

The power spectral density of the disc as a whole then becomes an incoherent sum of the 
PSDs of the individual disc radii, normalized by the square of the total luminosity of the disc. 
As an example, we take the instability time scales over which the flux varies to range from r = 
[a(/i/A) 2 f2] _1 to r = (af2) _1 , where as before a is the viscosity parameter, Q is the Keplerian 
frequency, and now the wavelength of the instability is given by A (h Ss A <C J"). (Note that for 
short wavelengths classical viscous and thermal instabilities merge into a single branch; Shakura & 
Sunyaev 1976.) As a simple example, let us take Ar ~ A oc Cr, where for our purposes C will be 
a dimensionless constant. For example, we might imagine that in a gas pressure dominated disk 
A ~ h ~ 0.02 yC 1 / 5 r- 21 / 20 ) giving C ~ 10~ 2 and nearly constant over a large range of radii. For this 
case we would be saying that the accretion rate is unsteady on radial scales less than a disc scale 
height, or equivalently the largest turbulent eddy size. We now write for the PSD of the disc: 



167r2 (a'fl)" 1 r 2 A 2 ?T 2 



, (27r//a'm 2 , 

w) — : . J — -2 . (5-4) 



(/ 



dr 2irr T 



where the sum is over radial bins of width Ar ~ A and a' is the constant of proportionality between 
the inverse instability time scale and the Keplerian frequency (equal to ~ a for thermal instabilities 
and/or short wavelengths and equal to ~ [h/X] 2 a for viscous instabilities). 

We analytically explore the behavior of this solution by substituting in the simple flux law 
T oc 7"~ 9 / 4 , along with the Keplerian frequency Q = c 3 /GM r~ z l 2 . Also, we replace the sum by a 
logarithmic integral over dr/Cr. This yields 

~ -j ; r - / t 1 23 dr > ( 5 - 5 ) 

4 a' c 6 J 6 1 + a J" 

where r is now in the units of GM/c 2 and a = (2-7T jGMj a! c 3 ). The above integral can be performed 
analytically; however, it is a bit more enlightening to look at the limits a <C 1 and d>l. In these 
limits, we obtain 

^ V2 Cj 2 GM _ 2/3 

b a' c* 5 . 

(5-6) 

1 Cj 2 GM 

a (a > 1) . 



576^6 a' c 3 
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The low frequency PSD has a slope that agrees well with observations (though we note that if we 
had chosen Ar oc C we would have found V(f) oc ln[l + a _2 /216] instead). If we set M = 6 Mq we 
then have 

C] 2 f f 



nf) 



0.15 



a 1 !/ 3 V 10- 2 Hz 



-2/3 



6 X 10~ b C j 2 a' 



f 



(a< 1) 



(a> 1) 



(5.7) 



10 2 H Z/ 

The low frequency solution has an amplitude in reasonable agreement with the observations for 
C ~ 0.01, a' ~ 0.1, and ~ 0.3. That is, flux variations on the order of 30% occuring over extents 
of roughly a scale height yield power comparable to the observations. Note, however, that with the 
same parameters the high frequency solution yields too little power. 
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Figure 2. Combination of viscous/thermal time scale and acoustic mode power compared to the 
observed PSD of GS 1124-683 (circles; cf. Miyamoto et al. 1993). In both graphs we have set 
M = 6 Mq, C = 0.1, a = 0.3, and a' = 0.03. As in Figure 1, the bandpass for both graphs 
was taken to be 1.2 — 15.7 keV. Left: Shakura & Sunyaev disk (cf. equations [4.21]-[4.28]) with 
Cj 2 = 1.3 Xl0~ 3 . Right: Flux is blackbody at each radius (cf. eq. [4.21]) with the peak temperature 
being artificially set to 1 keV and Cj 2 = 8 X 10~ 4 . 

In the above, we have used an oversimplified flux model and have not included a bandpass. 
As before, a bandpass in the soft X-ray range will exclude the outer regions of the disc, and hence 
lead to a reduction of the low frequency power. In addition, for accretion discs with "no-torque" 
inner edge boundary conditions the flux drops to zero on the inner edge of the disc. This will lead 
to a further reduction in the high frequency power. In Figure 2, we present numerical calculations 
wherein we use a bandpass of 1.2 — 15.7 keV and use the flux laws and disc models of §4. The 
left graph is for a Shakura- Sunyaev disc (cf. equations [4.21]-[4.28]) with M = 6 Mq, C = 0.1, 
and a = 0.3. We have also taken a' = 0.03 (which for a viscous instability would correspond to 
h/X ~ 0.3) and we have set Cj 2 = 1.3 X 10~ 3 (i.e. ~ 30% flux variations over radial extents the order 
of the disk thickness). The right graph is for the flux law of equation (4.21) (i.e. blackbody), with 
the peak temperature again being artificially set to 1 keV. Here also we have taken M = 6 Mq, 
C = 0.1, a = 0.3, and a' = 0.03. We have set Cj 2 = 8 X 10~ 4 . In each of these graphs we have 



16 



added in the dynamical power as presented in Figure 1 (for the same parameters). As before, the 
circles represent high state data of GS 1124-68, for which M ~ 6 Mq and C ~ 0.1. (Note that in 
this figure we have not attempted to find a rigorous "best fit" to the data. Our parameters were 
chosen by eye in order to give a general indication of the trends.) 

The combination of viscous and dynamical power is mildly successful in reproducing the data. 
The amplitude and rough shape of the low and high frequency data are reproduced, however, there 
are significant deviations. Note that the shape of the low frequency power for the Shakura-Sunyaev 
discs tends to be flat at very low frequencies, and then rolls over into the f~ 2 behavior. The expected 
/~ 2 / 3 behavior is absent. This is because the high frequency cut-off (due to the no-torque boundary 
condition) and the low-frequency cut-off (due to the bandpass) restricts us to a narrow range of 
radii. We are mostly seeing the V(f) oc 1/[1 + (27r/r) 2 ] behavior from this narrow range of radii. 
The blackbody flux law of equation (4.21) with the peak temperature set to 1 keV is observable 
over a slightly broader range of radii, so we see some of the expected /~ 2 / 3 behavior, which is in 
better agreement with the data. This model also has a broader high frequency peak, which also 
agrees better with the data. We see that in order to reproduce the data with viscous/thermal 
and dynamical time scales, we require a disc that emits soft X-rays over a broad range of radii. 
It is unlikely, however, that any realistic flux model will allow for much power above ~ 10 2 Hz. 
Furthermore, we expect that with most flux models the PSD must also flatten below ~ 10~ 2 Hz. 
Again, observations by XTE and USA will be able to confirm or refute this behavior. 

6 CONCLUSIONS 

If we are to look for observational clues that elucidate the basic physics of accretion discs we should 
explore the high state of galactic black hole candidates. The high state appears to be well explained 
by optically thick accretion disc models and appears to be relatively quiet and stable on short time 
scales. This is to be compared to the non-thermal and highly variable low state, which may require 
the addition of clouds, winds, and/or optically thin regimes in order to explain the observed energy 
spectra. In this work we have explored the consequences of the assumption that viscosity (and 
hence angular momentum transport and energy generation) is ultimately due to accretion disc 
turbulence. 

We explored two ways in which disc turbulence could modulate the observed emission. First, 
we considered the role of advection; however, we found that typical photon diffusion times through 
a scattering atmosphere are shorter than eddy rollover times. Thus eddies are incapable of "dredg- 
ing" hot photons from the disc midplane. Second, we considered direct modulation of the emission 
rate. Under the assumption that locally within the disc the energy generation rate scales as a small 
power of the density, one expects flux modulations to arise from the weak density fluctuations 
associated with the turbulence. These weak density fluctuations can be considered to be essentially 
acoustic modes driven by the turbulence. The amplitudes of these acoustic modes are proportional 
to the Mach number of the turbulence, which in our model is related to the familiar disc parameter 
a. Note that our model assumed isotropic hydrodyamic turbulence, whereas many recent models 
employ either two dimensional hydrodynamic turbulence or three dimensional magnetohydrody- 
namic turbulence (cf. the references of §1). For our model, at least the largest eddies (with the 
lowest frequencies) were capable of modulating the flux. 

We then considered how this variability would be manifested in the observed X-ray power 
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spectral density (PSD). We found that the theoretical PSD could be derived from summing the flux 
variations of individual modes over the mode phase space density. We compared these predictions to 
high state observations of GS 1124-683, which show the characteristic low amplitude V(f) oc f~ 07 
of BHC high states (cf. Miyamoto 1994). Our model was able to reproduce much of the observed 
high frequency power; however, when we included the X-ray bandpass of the detectors we were 
unable to reproduce the low frequency power. In addition, we required a large a (<; 0.3) in order to 
reproduce the amplitude of the high frequency power. This large an a is inconsistent with estimates 
of two dimensional hydrodynamic turbulence (cf. Narayan et al. 1994); however, it is as yet unclear 
whether or not this large an a is inconsistent with MHD turbulence as embodied in the magnetic 
shearing instability framework (cf. the references of §1). Note also that in equation (4.13) we have 
assumed that eddies of all sizes are capable of modulating the emission. We have not accounted for 
the fact that at some radii and depths only the largest eddies contribute to the modulation. Taking 
this detail into account would further reduce the high frequency power; however, this should not 
drastically effect the simple estimates presented in §4. 

Finally, we considered the possibility that the low frequency power could be caused by weak 
fluctuations on viscous/thermal time scales. We saw that flux variations on the order of 30% over 
distances comparable to the largest eddy sizes could account for the amplitude and rough shape of 
the low frequency power. However, viscous/thermal fluctuations were incapable of simultaneously 
reproducing the low and high frequency power. Combining the acoustic mode power with these 
crude fluctuation estimates was roughly capable of reproducing the observed PSD over three-four 
decades of frequency, though there were noticeable deviations. 

We note that it is difficult for any single variability mechanism to reproduce four decades of 
power if the variabilty time scale is proportional to the Keplerian time scale. The X-ray producing 
range of the disc exists between radii of r ~ 6 — 150 GM/c 2 . Since the Keplerian time scale 
goes as 7" 3 / 2 , any variability mechanism proportional to this time scale can cover only two orders 
of magnitude of frequency in the X-ray producing region of the disc. We see that two separate 
time scales can cover perhaps four decades, as we showed in §5; however, unless there is a third 
mechanism at work we expect there to be a break for / Ss 10~ 2 Hz. In addition, as the dynamical 
time scale is likely to be the fastest time scale in the disc we expect a sharp rollover for / <; 10 2 
Hz. For our model of turbulence, this rollover should be oc /~ 5 . We expect that both XTE and 
USA should be able to observe the low frequency break and high frequency rollover if our ideas are 
correct. 

We see two major areas in which the above model can be improved: the model of turbulence, 
and the model of emission from the disc. We expect that any significant progress made will come 
from attacking both of these problems simultaneously. Of great interest would be a more complete 
understanding (in a global sense) of the magnetic shearing instability, if this is indeed the correct 
explanation for accretion disc viscosity. One needs to understand the amplitude of the turbulence, 
how it couples to the energy generation rate, and the characteristic viscous time scale. Also, as 
we have seen from Figure 2, the particular flux law that we use affects the range of frequencies 
over which the the emission is modulated. It is unlikely, however, that any reasonable flux law will 
provide for dynamical variability power at frequencies much below a few Hz unless some mechanism 
such as Comptonization is invoked. We note though that the high state lacks an extended, hard tail 
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in the energy spectra, and therefore it is difficult to reconcile such mechanisms with the observations. 

Finally, we need to ask whether such a thing as accretion disk turbulence exists at all in these 
black hole candidate systems. To this end we note that the V(f) oc f~ 7 dependence appears to 
be a characteristic behavior of the thermal emission of BHC. Even when model fits indicate that 
less than 1% of the photon flux is emanating from a decidely non-thermal source, the V(f) oc f~ 7 
dependence remains - albeit at a very low amplitude (Miyamoto 1994). As the slope of the PSD 
implies, the integrated power must peak near the characteristic dynamical time scales of the disc. 
The exact nature of these characteristic oscillations are still in doubt; however, these weak variations 
seem to be fundamental to the emission mechanisms of this BHC state and perhaps are the best 
evidence for the existence of accretion disc turbulence. 
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